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Amplitude equations for 3D double-diffusive convection interacted with a horizontal 
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Three dimensional roll-type double-diffusive convection in a horizontally infinite layer of an un- 
compressible liquid is considered in the neighborhood of Hopf bifurcation points. A system of 
amplitude equations for the variations of convective rolls amplitude is derived by multiple-scaled 
method. An attention is paid to an interaction of convection and horizontal vortex. Different cases 
of the derived equations are discussed. 
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I. INTRODUCTION 

Physical systems, in which double-diffusion induced 
convection plays an essential role, are often found in na- 
ture. There are two components with significantly differ- 
ent diffusion coefficients in such systems. It can be heat 
and salt in the sea water, heat and helium in stellar atmo- 
spheres, or two reagents in chemical reactors. As a result 
of various spatial distribution of these components in a 
gravitational field arises convection, which can have vari- 
ous forms and lead to a variety of phenomena 1] . Widely 
known, for example, salt fingers are, arising in salted and 
warmed from above water. It is understandable that the 
results of double-diffusive convection, for example, in the 
ocean, can be applied to double-diffusive convection in 
astrophysical systems or in chemical reactor. 

There are a number of works devoted to various the- 
oretical models of systems with double-diffusive convec- 
tion. In 80-90 years the formation of structures in the 
neighborhood of Hopf bifurcation points for the horizon- 
tally translation-invariant systems was actively studied 
in some works. The development of oscillations in such 
systems give rise to different types of waves (eg, stand- 
ing, running, modulated, chaotic), which well described 
by a generalized Ginzburg-Landau equations [1, Q . The 
equations of this type must be derived from the basic sys- 
tem of partial differential equations for the given physi- 
cal system by asymptotic methods. However, a full and 
well-grounded derivation of amplitude equations for sys- 
tems with double-diffusive convection (especially three- 
dimensional) is still poorly represented in the literature. 

The purpose of this work is the derivation of amplitude 
equations for the three-dimensional double-diffusive sys- 
tem in the neighborhood of Hopf bifurcation points for 
the case of roll-type convection. This extends the idea of 
previous work [3, 4] , where a two-dimensional and three- 
dimensional convection in a square-cells was investigated 
by alike methods. 



II. FORMULATION OF THE PROBLEM AND 
BASIC EQUATIONS 



Consider 3D double-diffusive convection in a liquid 
layer of a width h, confined by two plane horizontal 
boundaries. The liquid layer is heated and salted from 
below. The governing equations in this case are hydro- 
dynamical e_2_uations for a liquid mixture in the gravita- 
tional field 



dtv + (vV)v = -p ^ Vp + vA-\ 
dtT + {^V)T ^ X^T , 
dtS+{wV)S = DAS, 
div V = . 



Where v{t,x,y, z) is the velocity field of liquid, 
T(i, X, y, z) is the temperature, S{t, x, y, z) is the salt con- 
centration, p{t,x,y,z) is the pressure, p(t,x,y,z) is the 
density of liquid, g is the acceleration of gravity, v is the 
kinematic viscosity of fluid, x is the thermal diffusivity 
of the liquid, D is the salt diffusivity. Cartesian frame 
with the horizontal x-axis and y-axis is used, while the 
z-axis is directed upward and t is the time variable. 

Distributed sources of heat and salt are absent. On the 
upper and lower boundaries of the layer the constant val- 
ues of temperature and salinity are supported, higher at 
the lower boundary. The governing equations are trans- 
formed into dimensionless form with the use of Boussi- 
nesq approximation and following units for length, time, 
velocity, pressure, temperature and salinity respectively: 
h, h'^/Xy x/h, pox^/h?, Ta, S'a, where Ta and 5a are 
temperature and salinity differences across the layer. 
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The dimensionless governing equations for momentum 
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and diffusion of temperature and salt are [J] : 

ut + {uux + vuy + wuz) ~ —px + uAu, 
vt + {uvx + vvy + wvz) = -py + aAv, 

Wt + {uWx + VWy + WWz) = 

- Pz + aAw + (7Rt9 - aRsS,, (1) 
Ot + {u0x + vOy + w9z) -w^A9, 
6 + {u£,x + viy + w4) - w = rA^, 

Ux + Vy + Wz = 0. 

Where a = vjx'is the Prandtl number (a w 7.0), r = 
D/x is the Lewis number (0 < r < 1, usuaUy r = 0.01 — 
0.1). Rt — {get I xv)T i\ is the temperature Rayleigh 
number and Rs = {g^' h'^ / x^) S a is the sahnity Rayleigh 
number. 

Free-slip boundary conditions are used for the depen- 
dent variables (the horizontal velocity component is un- 
defined): 

Uz — Vz~w = 6 = £^~QdXz — 0, 1. 

It is believed that they are suitable to describe the con- 
vection in the inner layers of liquid and do not change 
significantly the convective instability occurrence crite- 
ria for the investigated class of systems 0. 



III. DERIVATION OF AMPLITUDE 
EQUATIONS - GENERAL FRAME OF 
DECOMPOSITION 

Consider the equations for double-diffusive convection 
in the vicinity of a bifurcation point, the temperature 
and salinity Rayleigh numbers for which are designated 
as and R*g respectively. In this case the Rayleigh 
numbers can be represented as follows: 

i?T-i?T(l + e'r-T), Rs^R*s{l + e''rs). 

Values of tt and rg are of unit order, and the small 
parameter e shows how far from the bifurcation point 
the system is. To derive the amplitude equations we use 
the derivative-expansion method, which is a variant of 
multiple-scale method [1, [§]. Introduce the slow vari- 
ables: 

Ti = et, Ta = eH, Xi = ex, Yi = ^y. 

In accordance with the method chosen, we assume that 
the dependent variables now depend on t, Ti, T2, a;, y, z, 
Xi, Yi, which are considered independent. Also replace 
the derivatives in equations ([1]) for the prolonged ones by 
the rules: 

dt^ dt + edxi + e^^Ta , dx ^ dx+ edxi , dy \/e% . 



Then the equations ([T]) can be written as: 

Ut + {uUx + WUz) +Px - iy{Uxx + Uzz) = --/evuY^ 

- e[uTi + uux^ + pxi - "^aUxXi - cuy^yj 
-I- s^[auxiXi -UT2], 

Vt + [uVx + WVz) - a{Vxx + Vzz) = -V^ivVYi + PYi) 

- e[vTi + uvxi - 2aVxXi - o-vy^Yi] 
+ e^[avxiXi -VT2], 

Wt + {uWx + WWz) +Pz - 0-{Wxx + Wzz) 

- gRtO -\- (tRs^ = —\/evwYi 

- e[wTi + uwxi - 2aWxXi - a-wy^Yi] 
+ e^[crwxiXi - WT2 + crrrO - crrsS,], 

Ot + [uOx + wez)-w~ {Oxx + Ozz)^ -VIvOy^ 

- e[0T^ + uOx, ~ 2exX^ - ^Y^Y^] 
+ ^^[OxiXi - ^Ta], 

6 + {uix + wS,z) - w - t{S,xx + izz) = -\/ev^Yt 

- e^Ti + <Xi - 2t^xX^ - rCi-iYi] 

Ux + Wz^ -\/evYt - ewjfi • 

(2) 

We seek solutions of these equations in the form of 
asymptotic series in powers of small parameter e: 

u = eui + e^U2 + £^uj, + • • • , 

V = e^/evi + e^\fev2 + e^Vevy, -)-•••, 

w = ewi + e'^w2 + e^wa + • • ■ , 

9 (o) 
p^epi+e P2 + £ P3 H , 

e = e9i+ 6^02 + 6^63 + • • • , 
e - eCi + + + ■ ■ ■ ■ 

After their substitution in ^ and collection the terms at 
e" we obtain the systems of equations to determine the 
terms of the series ([3]) . 

IV. TERMS OF THE FIRST ORDER IN s 

At O(e^) we obtain the following system: 

Ult + Plx - Cr{uixx + Uizz) = 0, 
Wit + PlYi - 0-{vixx + Vizz) = 0, 

Wit + Piz - o-{wixx + wizz) - (jRtOi + a-Rs^i = 0, 

Olt-Wi-{0ixx+9lzz)=O, ^ ' 

Clt -Wi- t{S,Ixx + Clzz) = 0, 
Ulx + Wiz = 0. 
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Solutions in the form of normal modes (roll-type cells) 
are: 



Wi 



Ae 



""^ e^' sinTTz + cc, 



h = Ola e"""^ e^* sinTTz + c.c. 



""^ e^' sinTTz + c.c, 



Pi 



Pia e''""' e^* cos TTz + pi 



ui — uia e''^^ e"^* cos irz + iti 



c.c, 
- c.c. 



(5) 



Via e'*^"^ e"^' cos ttz + vi + c.c. . 



Terms with the bars and hats depend only on the slow 
variables Xi,Yi,Ti,T2. Substituting the expressions ([5]) 
in the equation ^ gives relations for the amplitudes. 
Hereinafter = fc^ + tt^ is the full wave number. 
Parameters A, k, Rt, Rs are related by: 

+ a{k^/>c^)[Rs{\ + >c^)-RT{X + TM:'^)] =0. 
This equation has three roots, two of which can be com- 
ics I 




II 



Rt 



Here L is linear differential operator such that Lipi = 
corresponds to the system Qj, where — 
{ui,Vi,Wi,Oi,^i). Functions Qi include terms, resonat- 
ing with the left side of the equations, namely: Qi = 
^ + Qf^ + Qf^- Here Q^p and Qf^ generate secu- 
lar terms of the two types in solutions, and Q\ does 
not generate secular terms of any kind. The condition of 
absence of secular terms of the first type consists in the 
requirement of orthogonality of functions Q^f'^ and the 
solution of the adjoint homogeneous equation L*ip* — 
, which usually takes the form of an amplitude equa- 

tion. Terms Qi are constants with respect to the fast 
variables, and to exclude the violation of the regularity 
of the expansions ([5]), they should be equated to zero 
These conditions also have the form of amplitude equa- 
tions. 

As a result, the condition of absence of secular terms 
of the first type in solution of equations at is as follows 
(q;o is defined in ([7]) below): 



Ati + ao{2ikAx^ + ^nyj + ikAVL 



To exclude secular terms of the second type, we introduce 
the horizontal stream function f2 such that: 



■^1 



VI. ^n-SYSTEM OF AMPLITUDE EQUATIONS 



FIG. 1: Plane of the Rayleigh numbers Rt — Rs- On the 
line 1 takes place Taylor bifurcation, on the line 2 takes place 
Hopf bifurcation. 



We write the resulting amplitude equations for the sys- 
tem at e^: 



plex conjugates. In the case of Hopf bifurcation these 
two roots acquire positive real part at some R^ (to is a 
frequency of convective waves): 



JXt 



C + T „^ >i , 

-R^ + —j{l + r){r + CT), 



1 - T 
1 + 



(jR*s^-t''k\ 



V. THE EQUATIONS AT AND 



At2 =rA + ai \^dxi + ~ aoAx^Xi 

+a2A\A\'' + J{fl,A)+F{fl,A), 
{nr, - (T^x^x, )x, - J{n, nx, ) + GiA) . 



(6) 



Here the Jacobian J(f7, /) = {Qxi/yi — ^Yifxi) is intro- 
duced. Also operators A) and G(A) are defined as 
follows: 



The obtained systems can be written in general form 



as: 



Lipi = Qi 



F(n,A) =a3{ikAnx,Y, + Ay.nY.Y,) + a^Any.Y.Y,, 
G{A) = pd^nHn + ^Re{ik{AA*yJx,). 
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FIG. 2: Rc(a3(w)) (left) and Im(Q3(w)) (right) with different FIG. 3: Re(a4(w)) (left) and Im(a4(w)) (right) with different 
k: k = 1.75 (dots), k = tt/v^ (solid), k = 3 (dash). At cr = 7 k: k = 1.75 (dots), k = tt/n/S (solid), k = 3 (dash). At cr = 7 
and T = 0.02. and t = 0.02. 



The coefficients of these equations are given by the ex- 
pressions: 



1 ( 



(1 + T + a)x- 
11(1- 



Q!2 



a4 



4ia; ' 

ifc 



.2 



(7) 



Q!3 



2ifc 



1+1 fc2 -l)/?2-/3/33 



(1 - /32 - /3/33) . 



Here for the convenience and compactness of the formulas 
we introduced functions: 



/3 = 

/32 
/?3 



2fe2 



K (r + (T + TiT)ia; + rcrx 

~~2 ' 



iw + (1 + T + C7)x2 



(r + cr + TC7)ia; + rax^ 
iw + (1 + T + cr)x2 ' 

2ia;(iw + (1 + r + cr)^^) ' 

2iw(iw + x2)(ia; + rx^) 

x^((l + T + 2cr)ia; + (1 + + tct + (7)x2) 
2(ia; + x2)(ia; + Tx2)(ia; + (1 + r + C7)x2) 



It is easy to see that for k = 7r/\/2, which corresponds 

to the fastest growing mode for not too high frequencies, 
the parameter /3 vanishes, and the above formulas are 
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significantly simplified. 



^ ~ \6iuj) 2iw + 3(l + T + cr)7r2 



ao 



as 



VII. 



2iw 
i2V2 



ai 



At: \uj 



Pi, 
i%/2 



a2 



8iw 



(l + /?2), a4=^ (1-/32) 



COMPLEX GINZBURG-LANDAU (CGL) 
EQUATION 



Let us consider some special cases of the system ([6]) to 
which it reduces at various simplifying assumptions. The 
simplest form of the equation is obtained if we assume 
that we are considering: 1). neglect the interaction with 
the vortex, and 2). consider the dynamics on the single 
spatial variable, and 3). believe that the wave number 
corresponds to the first losing stability mode k — 
Under these conditions, the system takes the following 
form: 



At^ =rA + asAx.x, + a2A\Af 



(8) 



Coefficients of the obtained CGL equation (consistent 
with are: 



^5 = ^ + -■ — Pi = — + -■ — 



(r + (7 + rcr)ia; + rayp' 
icj + (1 + r + cr)>f2 



In the case when the wave number k is not restricted 
by any conditions the coefficient as is expressed in gen- 
eral terms as as = ai — ao (see fig. S]). In the limit of 
high frequencies lo the resulting equation reduces to the 
nonlinear Schrodinger equation and has such solutions as 
"dark" solitons [3]- 



VIII. CONCLUSION 

• The Afi-system of amplitude equations ([6]) describ- 
ing 3D double-diffusive roll-type convection inter- 
acting with horizontal vorticity field was derived. 

• An approach to calculation of amplitude equation 
coefficients which allows to get relatively compact 
formulas such as ^ was developed. 

• As a special case the complex Cinzburg-Landau 
equation ([5]) for the 2D double-diffusive convection 
for an arbitrary k (width of the convcctivc cells) 
was derived. 




15 20 25 30 



FIG. 4: Re(ci5(tj)) (left) and Im(a5(i.^)) (right) with different 
fc: fc = 1.75 (dots), k = 7r/\/2 (solid), it = 3 (dash). At cr = 7 
and r = 0.02. 
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The results can be used to describe the processes of heat 
and mass transfer, the formation of vortex structures in 
the ocean and the atmosphere by convection, and may 
also be the basis for constructing more advanced models 
of this kind. 
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